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Monojet events at colliders have been used to probe models of dark matter and extra dimensions. 
We point out that these events also probe extensions of the Standard Model modifying neutrino- 
quark interactions. Such nonstandard interactions (NSI) have been discussed in connection with 
neutrino oscillation experiments. Assuming first that NSI remain contact at LHC energies, we 
derive stringent bounds that approach the levels suggested by the 8 B solar data. We next explore 
the possibility that the mediators of the NSI can be produced at colliders. The constraints are 
found to be strongest for mediator masses in the 10 2 - 10 3 GeV range, with the best bounds above 
~ 200 GeV coming from ATLAS and below from CDF. For mediators with masses below 30 GeV 
the monojet bounds are weaker than in the contact limit. These results also directly apply to light 
dark matter searches. Lastly, we discuss how neutrino NSI can be distinguished from dark matter 
or Kaluza-Klein states with charged lepton searches. 



I. INTRODUCTION 

Many extensions of the Standard Model predict new 
weakly interacting particles that, if produced at colliders, 
would escape the detector leaving an apparent imbalance 
of energy and momentum. Particularly striking events of 
this type are dominated by a single energetic jet recoiling 
against "nothing". Well-studied searches utilizing these 
so-called monojet events include the ADD-type models 
of extra dimensions, in which the invisible new physics 
("nothing") is Kaluza-Klein gravitons [IH5], and models 
of dark matter, in which the latter is produced directly 
from colliding partons [tjHTU]. 

The main "irreducible" Standard Model (SM) back- 
ground to the monojet searches - especially at high trans- 
verse momenta of the jet - is provided by neutrinos, 
which are created in the decays of the Z bosons, or the 
W bosons (when the accompanying charged leptons are 
missed). Yet, neutrinos are not necessarily just a nui- 
sance to new physics searches. As we discuss in the 
present Letter, neutrinos themselves could be affected 
by new physics modifying their production rates. "Non- 
standard" interactions (NSI) of neutrinos could thus fake 
the signal of dark matter or extra-dimensional physics. 

The idea of neutrino NSI is, in fact, not new. It is 
prominently featured already in the seminal paper by L. 
Wolfenstein [TT] , which laid the foundation for the MSW 
effect [TTJ H2] ■ Hundreds of subsequent papers explored 
the oscillation impact of NSI in various scenarios. Dur- 
ing the 1980's and early 1990's, due to limited available 
data, NSI were mainly discussed as an alternative mech- 
anism to the mass-induced oscillations (e.g., Q21 HI])- 
This changed in the last decade, thanks to the dramatic 
advances in solar, atmospheric, reactor and beam neu- 
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trino experiments. It is now possible to search for rel- 
atively small, subdominant effects in oscillations caused 
by NSI (see, e.g., |15j for an overview). 

To illustrate this, we plot in Fig. [I] the survival prob- 
ability P(y e — > v e ) of solar neutrinos, with only SM 
physics (thick black curves) and with the addition of 
NSI. Here, the neutrino-quark NSI couplings e qP (see 
our conventions in Eq. ([l]) below) are a few percent of 
the SM weak interactions. We see that the shape of the 
survival probability at E v ~ a few MeV, in the transi- 
tion window between vacuum (low energies) and MSW 
(high energies) regimes, is a sensitive probe of neutrino- 
matter interactions. This is especially so when the flavor- 
changing component is introduced (right panel), as has 
been noted before (c/. Fig. 1 in [18]). 

Curiously, the latest solar neutrino results are better fit 
with NSI than with the SM interactions alone. The SNO 
[2~Tl |2"2"] and Super-Kamiokande [22] experiments both 
lowered their energy thresholds in recent years, aiming 
to observe the standard MSW "upturn" of P(v e — > v e ) 
in the transition window. Yet, neither experiment has 
seen it. Additionally, the Borexino experiment targeted 
8 B neutrinos [24] in the same energy window and likewise 
found no upturn. A careful analysis of the combined data 
a year ago [21] found that nonzero NSI were favored at 
the ~ 2ct level, a result that may strengthen with the 
addition of the recent SNO data [22] . 

Can other data exclude the NSI couplings in the range 
favored by the solar 8 B data? It turns out that the Super- 
Kamiokande atmospheric neutrino data [19] do not, even 
with the addition of the recent data from MINOS |20j . 
Neither do a variety of other experiments that are sensi- 
tive to neutrino NSI [14j [26H28] , at least not in a model- 
independent way. 

Could the LHC and Tevatron monojet datasets be 
more sensitive to the neutrino NSI that the solar neu- 
trino data? As we show in what follows, the answer de- 
pends on the scale of new physics. We present bounds 
for different assumptions about this scale. 

As already mentioned, the monojet signatures of neu- 
trino NSI may look exactly like those of light dark mat- 
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FIG. 1: The effect of the flavor-diagonal (left) and flavor off-diagonal (right) NSI on the day-time survival probability P(y e — > v e ) 
of electron neutrinos from the Sun. The thick black curves represent the Standard Model expectations, using the recently 
measured sin 2 26\s ~ 0.1 [161 117] - while the thin red curves represent the result of varying the NSI e parameters per electron 
in the range [—0.2,0.2]. The neutrino is taken to be produced at the center of the Sun (a good approximation for the 8 B 
neutrinos). 



ter or large extra dimensions. In fact, our monojet con- 
straints can be directly recast as bounds on dark matter 
scenarios. With the addition of other data, however, it 
may be possible to resolve the "dark matter/neutrino 
ambiguity" . We present several examples of this. 

The presentation is organized as follows. After a brief 
introduction (Section [Tip , we analyze the potential of 
monojet searches under the assumption that neutrino 
NSI remain contact (Section III). We show that the 
present data allow the scale of these operators to be as 
low as 500 GeV, which motivates us to consider scenar- 
ios with finite mediator mass in Section |TV} In Section [V] 
we discuss how multilepton searches at the LHC as well 
as lepton flavor-violating decays can be used to discrimi- 
nate neutrinos from other sources of missing energy. Sec- 
tion IVl] summarizes our conclusions. 



II. GENERALITIES AND NOTATION 

We begin by defining the Lagrangian for neutrino 
NSI. We consider modifications to the neutral current 
neutrino-quark interactions. The strength of these modi- 
fications is conventionally defined in units of the SM weak 
interaction, given by Gp- 
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Here / denotes the SM fermion flavor, P is the left/right 
projector, and e^ p are hermitian matrices in the neutrino 
flavor space spanned by a,/3 = e,/z,r. Throughout the 
Letter we assume that the neutrinos are left-handed and 
consider / = u, d. The up and down quark couplings 
are relevant for neutrino oscillations in matter and also 
provide the dominant contribution to proton collisions. 

Importantly, the relationship between NSI effects in os- 
cillations and at colliders is not one-to-one. Indeed, since 



for oscillations in matter forward scattering amplitudes 
add up coherently, only vector couplings + are 
important. In contrast, in collisions nonstandard axial 
couplings also modify the neutrino production rate and 
hence are also probed. Moreover, note that NSI in Fig. [TJ 
and in many oscillation analyses are given per electron. 
Since, for the chemical composition of the Sun, there are 
4-5 quarks per electron, the range of the NSI parameters 
in Fig. [T] is a few percent per quark. 

As the right panel of the figure shows, flavor-changing 
NSI of this magnitude (and the right sign) make the 
P(v e — > v e ) function above a few MeV flat [TH[. This 
fits the data from SNO [23 H2|, Super-Kamiokande [23], 
and Borexino [23] better than the SM curve [25] . 

The implicit assumption in Eq. ([IJ is that the new 
physics can be safely integrated out, leaving a contact in- 
teraction. This seems reasonable at energy scales relevant 
to solar neutrinos. In the neutrino oscillation literature, 
this assumption is also typically extended to the more 
energetic atmospheric neutrinos, where it is less obvious. 
At the Tevatron and LHC energies, it becomes even less 
obvious. We will therefore explore the collider signatures 
of NSI in two stages: first, by assuming the contact form 
of Eq. (IT]) and then by relaxing this assumption. 

Eq. (TTh in general contains both flavor-changing and 
flavor-diagonal NSI. The former produce final states that 
have no SM analogs, and hence behave at colliders like 
light dark matter. In contrast, the latter can interfere 
with the SM, leading to a nontrivial difference with the 
dark matter analyses. Whether this interference is practi- 
cally important dependends on the strength of the bound, 
as we will explore in what follows. 

Another important difference with dark matter is that 
neutrinos are charged under the electroweak symmetry. 
This suggests that NSI may be accompanied by same 
strength operators involving the charged leptons. This 
is indeed so if before electroweak symmetry breaking the 
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interactions leading to ([I]) can be written as the following 
dimension-6 operators 



/»dim- 
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2e qP _ 

—^(L a rLp)(q^Pq), (2) 



where L = [y, t) is the lepton doublet and v 2 — 1/ \/2Gf- 
These operators are very strongly bounded by processes 
involving charged leptons £. It has been argued, how- 
ever, that Eq. ([2| should not be used to derive model- 
independent bounds, as the NSI could also arise from 
more complicated effective operators. If such operators 
involve the Higgs field, the obvious SU(2)l connection 
may be broken [H] [2&H28] . Typical examples are mod- 
els where arises from dimension-8 operators of the 
form 27J 
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TABLE I: Bounds on the contact NSI from the CDF and 
ATLAS monojet + MET searches. The CDF bounds are 
based on 1.1 fb _1 of data and are shown for two sets of cuts, 
the softer "Generic Search for New Physics" (GSNP) cuts [32] 
and the harder ones optimized for the ADD searches [HE]. 
The ATLAS bounds are based on 1 fb -1 for the three different 
cuts analyzed in [31] . All bounds correspond to 95% C.L. The 
bounds do not depend on the neutrino flavor a, /3 = e, fi,r nor 
on the chirality P — L, R of the quark. We assume only one 
coefficient at a time is turned on. When several coefficients 
contribute the bound reads as shown in Eq. 



with H being the Higgs doublet. In defining the coeffi- 
cient of the operator we used the fact that in the unitary 
gauge W H —> (v + h) 2 /2, with h the Higgs field. In 
this case the low-energy Lagrangian ([I]) need not be ac- 
companied by same-strength operators involving charged 
leptons. 

Lastly, let us note that even the NSI Lagrangian ([3| 
will inevitably contribute to charged lepton processes at 
high energies 129] . We will see in Sec. VB that the op- 
erator in Eq. Q does indeed produce charged leptons at 
the LHC, at potentially detectable levels. 
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FIG. 2: Feynman diagrams contributing to the monojet sig- 
nal Q, with time flowing from left to right. The shaded blobs 
denote the NSI contact interaction. At the 7 TeV LHC the qq 
initial state contributes approximately the 70% of the signal. 



III. MONOJET BOUNDS ON NEUTRINO 
CONTACT INTERACTIONS 

At the simplest level, the four fermion operator in 
Eq. gives rise to the distinctive but invisible pro- 
cess qq — > v a vp. This event is rendered visible if for 
example one of the initial state quarks radiates a gluon, 
qq — » V a vpg. This along with the two other diagrams in- 
volving quark-gluon initial states shown in Fig. [2] consti- 
tute the monojet plus missing transverse energy (MET) 
signal we consider here: 

PP ipP) -> 3 v a vp , j = q, q, g. (4) 

Analogous constraints on NSI [27] and dark matter [5U] 
involving electrons arise at e + e _ colliders where instead 
of a jet one has a photon in the final state. 

Below, in Sec. |III A[ we describe our derivation of 
the bounds from the LHC (ATLAS [3T]) and Tevatron 
(CDF [4, 5, 32 ) data, assuming the interactions remain 
contact for all relevant energies. The summary of these 
bounds is presented in Table [T] We note that these con- 
straints improve considerably the corresponding bounds 
on £ eT , e TT , e ee , as reported in [28] , 

Given that the LHC is already at the frontier of 
neutrino-quark interactions, it is natural to ask how these 



bounds will change in the near future, as more data is 
collected and analyzed. In Section |IIIB| we attempt to 
make some informed projections of the bounds, conclud- 
ing that a significant improvement in the bounds will only 
be achieved once systematics are reduced. We note that 
although CMS also has a monojet study with a compa- 
rable data set [33], we use the ATLAS study precisely 
because of its careful discussion of the systematics. 

We also examine the effect of the event selection crite- 
ria as a determinant in setting the bounds. In particular, 
note that while the hardest px cut of the five selection cri- 
teria in Table [I] yields the strongest bound in the contact 
limit, the same is not true in the light mediator regime, 
as we show in Sec. IIVI 



A. Analysis details 

The standard model (SM) monojet backgrounds are 
primarily due to pp(pp) — » jZ — > jvv, pp(pp) — > jW -> 
j£v where the charged lepton is missed, and multi-jet 
QCD events [3TH33]. 

The CDF collaboration released its monojet data with 
two sets of cuts. One is designed for a generic search for 
new physics (henceforth, the GSNP cut) [32], the other 
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is specifically optimized for ADD searches [U [5] (hence- 
forth, the ADD cut). In the first case, the cut on the 
transverse momentum of the leading jet is rather mod- 
est, pt > 80 GeV; the missing energy is required to be 
> 80 GeV and the transverse momenta of the second and 
third jets (if any) have to be below 30 GeV and 20 GeV. 
In the second case, the cut on the transverse momentum 
of the leading jet is harder, pr > 150 GeV; the missing 
energy is required to be > 120 GeV and the transverse 
momenta of the second and third jets have to be below 
60 GeV and 20 GeV. 

ATLAS considered three different selection criteria re- 
ferred to as LowPt, HighPt, and veryHighPt cuts. The 
main difference between these is the cut on the transverse 
momentum of the leading jet, that respectively reads 
p T > 120, 250, 350 GeV. We also imposed the additional 
jet vetoes and further cuts as described in [31]. The to- 
tal systematic and statistical uncertainty amounts to ap- 
proximately 5%, 7%, and 13% of the predicted events for 
the three cuts considered. In addition, the uncertainty is 
dominated by systematics, as we discuss in some detail 
below (Sec. |nTB |. 

We generated the parton-level signal Q for a given set 

a, (3, f, P with Madgraph/Madevent_v5 [34]. The relevant 

Feynman diagrams for monojets from NSI are depicted 

in FIG. [2] We imposed a 50 GeV generator-level pt cut, 

and then passed the data to Pythia 8 [35] for initial and 

final state radiation, hadronization, and event selection 

and to Fast jet 2.4.4 36J for jet clustering. Multiple 

interactions were switched on and off and found not to 

affect our results. We have also explicitly checked that we 

do not double-count jets. By generating the parton-level 

process pp{pp) — > V a vp and allowing Pythia to ge nera te 

the jet, we find consistent results (here and in Sec. IV). 

f p — 

An upper bound on the coefficient e a p is found by 
requiring that the number of events that pass the cuts 
be below the 95% CL bound reported by the collabora- 
tions. From Table |T] we see that the LHC has already 
superseded the Tevatron in sensitivity to contact NSI. 
We further note that the ADD-optimized cuts used by 
CDF turn out to be suboptimal for the NSI search. 

As noted above, unlike dark matter monojet searches, 
flavor-diagonal NSI interfere with the SM. Turning on 
only the cross section for MJ) can be written as 
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Interference plays a significant role only for sufficiently 
small e Qa 's. For the bounds given in Table [I] we find 
interference to be subleading, implying a correction of 
less than ~ 10% to our bounds. For example, for the LHC 
at 7 TeV the up- type quarks give cr^fj = 1.2 pb, while 
interference contributes af^ — 2.6 x 10~ 2 pb, af^ = 
-5.9 x 10~ 2 pb. 

For off-diagonal couplings, note that once one of the 
e;^ is turned on the NSI operators generate not only (JiJ 
but also its conjugate pp — » jvpv a . These processes in- 
coherently contribute to the j+MET signal. Hence, the 
cross section a(pp — > j+MET) is effectively enhanced by 



a factor of 2 compared to the case of diagonal couplings. 
This leads to an improvement of a factor of v2 of the 
bounds, as shown in the last two lines of Table [I] 

Furthermore, though the bounds do not depend on the 
chirality P = L, R of the incoming parton, they are sensi- 
tive to the quark flavors / = u, d of the operators via 
the parton distribution functions. At both the LHC and 
the Tevatron the processes involving up-type quarks are 
enhanced, and the bounds on e uP are therefore stronger 
than those on e dp . 

Finally, we emphasize that the constraints reported in 
Table [T] apply when only one NSI coefficient is switched 
on at a time. More generally, however, the bounds can 
be summarized as: 
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Here, the flavor off-diagonal e's are to be summed twice, 
as in \ei P \ 2 + \e d T P \ 2 = 2\e dp \ 2 . The interference effects 
have been neglected, for the reasons explained above. 



B. Systematic Uncertainties and Projections 

An inspection of ATLAS's 31j Table 1 reveals that the 
dominant source of uncertainty for monojet searches at 
the LHC is due to systematics. Although most of this 
uncertainty (including jet energy resolution, parton dis- 
tribution functions, etc.) will presumably improve with 
statistics, it is clear that a luminosity upgrade will not 
lead to a simple y/N rescaling of the bounds. 

It is indeed precisely the dominance of systematic er- 
rors that make ATLAS's hardest p? selection better 
suited to constraining NSI contact interactions. In the 
absence of systematic errors, a % 2 statistic formed out 
of the signal and dominant Z — > vV background peaks 
at lower p^, implying that softer momentum cuts pro- 
vide more stringent bounds. When systematics are intro- 
duced, however, the significance of the signal is always 
reduced compared to the idealized statistics only case, 
and the optimal bound is obtained at the veryHighPt se- 
lection cut. In the absence of detailed knowledge of how 
the systematics vary with pt it is impossible to know if 
an even harder cut on the transverse momentum of the 
jet would lead to even more stringent bounds. 

Thus although we cannot obtain quantitatively pre- 
cise projections, it is clear qualitatively that the bounds 
will not change appreciably with luminosity unless the 
systematic errors are reduced. For example, using the 
X 2 statistic again, we find that even with 15 fb _1 at 
the 7 TeV LHC and with a factor of 3 improvement 
in the systematic uncertainty, the epsilon bounds of Ta- 
ble [I] are improved by less than a factor of 2. We there- 
fore conclude that the bounds in Table U will remain the 
strongest bounds for contact neutrino-quark interactions 
until a considerable reduction of systematic uncertainties 
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FIG. 3: Feynman diagrams contributing to monojet signals in 
the leptoquark model. Dashed lines denote a leptoquark. The 
last diagram can dominate for light leptoquark masses, but is 
subdominant at low energy as it leads to a dimension-8 op- 
erator involving a gluon field, two quarks, and two neutrinos. 
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FIG. 4: Monojet signals in the Z' model of the NSI contact 
operator. Wavy lines denote a Z' . 



to 




M z [GeV] 

FIG. 5: Contours of fixed generator-level cross-section in the 
Z 1 model. Here it is assumed that the Z' couples equally to 
ul and a flavor non-conserving neutrino pair. The red-dashed 
curve illustrates the naive bound obtained by using a fixed 
acceptance, corresponding to the contact-operator with very- 
High Pt cuts. See text for additional details. Actual bounds 
are shown in Fig. [7] 



is achieved . 

Finally, using the same v 2 procedure we can obtain 
a rough estimate of the bounds expected from the 14 
TeV LHC in an optimistic and completely unrealistic sce- 
nario where systematics are negligible. With 100 fb _1 of 
data at the 14 TeV LHC the bounds can be as strong as 
el P / P < 10" 3 . 



IV. MODEL-DEPENDENT BOUNDS 

The effective operator analysis of the previous Section 
presupposes that the scale of new physics is much higher 
than the energies probed in collisions. What happens 
when this is not the case? In this Section, we examine 
a scenario with a finite-mass mediator. We show that 
the contact limit does not set in at the LHC unless the 
mediator mass is above several TeV. We also find that, for 
very light mediators, the NSI parameters e are actually 
less constrained by monojets than in the contact limit. 



1 Note added: This point is illustrated perfectly by the recently 
released CMS monojet analysis with 4.7 fb — 1 of data I39| . This 
analysis employed event selection criteria similar to the very- 
HighPt ATLAS analysis |31| and arrives at similar systematic 
effects. The bounds we obtain from this CMS data are indeed 
essentially the same as the veryHighPt ATLAS bound shown in 
Fig. 7, despite five times more statistics. 



Any discussion beyond the effective operator limit is by 
necessity model-dependent. The effective operator of the 
form given in Eq. ([3| could be UV-completed in different 
ways. As an example, consider a i-channel completion 
with a leptoquark exchange between a quark and a neu- 
trino. The leptoquark in question is for example an elec- 
troweak doublet, color triplet scalar S with hypercharge 
Y = 1/6 that couples to the SM fermions via oIrLS. 
Higgs VEV insertions on the leptoquark line can account 
for a suppression of charged lepton processes The 
leptoquark would contribute to the monojet production 
rate via the diagrams shown in Fig. [3j It is instructive 
to consider how the NSI parameters in this model can be 
constrained by the wealth of the available data (beyond 
monojets). We will discuss this later (in Sect. VA). 



As a second example, consider an s-channel UV com- 
pletion with a Z' intermediate state. It is assumed that 
in the full model the SU(2)l symmetry is again appro- 
priately broken by Higgs VEV insertions. How exactly 
this is realized will dictate what other searches could be 
used to probe this scenario. For our immediate purpose, 
we are interested in the direct monojet bounds and hence 
will consider a schematic Z'-neutrino and Z'-quarks cou- 
plings (for an existence proof see |46| for an explicit 
model where the Z' couples only to quarks and neutri- 
nos but not to charged leptons). The relevant processes 
are shown in Fig. [2] These examples illustrate potential 
connections between neutrino NSI and various ongoing 
searches at the LHC. 
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Other models could be given (see, e.g., [37J[38]). Our 
goal here, however, is not to survey multiple specific sce- 
narios of new physics, but simply to demonstrate that 
the monojet bounds on NSI could vary significantly as 
a function of the mediator mass. To this end, we will 
specialize to the Z' model, and show how the monojet 
rates depend on Mz> and the coupling gz>- 

To begin, we compute the parton-level cross sections 
of the monojet process as a function of Mz> and gz' ■ 
For simplicity, the width of the Z' is calculated here as- 
suming coupling only to one quark flavor and chirality 
as well as one neutrino flavor, Y z i — gz'Mz 1 /67T ■. Wc 
consider proton-proton collisions at 7 TeV, and also spe- 
cialize to a flavor-changing NSI, so that the interference 
effects are absent. We again use Madgraph/Madevent_v5, 
which we set up to loop over a two-dimensional logarith- 
mically spaced grid of points. 

The resulting contours of constant parton-level cross 
section are shown in Fig. [5] The results are presented in 
terms of e = E^p^a (cf. Table|lj). We see here four regimes 
of interest: (1) the heavy mass, small coupling regime 
where the e cross section are independent of the media- 
tor mass, thus merging with the contact operator results; 
(2) the heavy mass, strong coupling regime (shaded tri- 
angle); (3) the intermediate mass regime where the cross 
section for fixed s is maximal; and (4) the low mass 
regime where for fixed e the cross section decreases as 
the mass is lowered. In all, we see that for a fixed value 
of the e parameter (fixed effect in neutrino oscillations) 
the monojet cross sections are indeed strongly sensitive 
to the mass scale of the mediator, varying by several or- 
ders of magnitude in the mass range [1, 10 5 ] GeV. 

The first regime (high mass, small coupling) is self- 
evident. In the second regime, the coupling gz> = 
y/2e(Mz' /v) becomes strong, the Z 1 becomes a very 
broad resonance, and the tree-level MadGraph treatment 
is clearly inadequate. In the third (intermediate mass) 
regime the mediator mass Mz' is of the order of the 
parton-parton collision energy. Monojet processes oc- 
curring via s-channel exchange are resonantly enhanced, 
compared to the contact regime. 

Lastly, consider the fourth regime, in which the cross 
section decreases as the mediator mass is lowered. In 
Fig. [5] this occurs for masses below a few hundred GeV. 
This happens because the typical momentum transfer 
flowing into the propagator dominates over the media- 
tor mass. In this limit for fixed coupling gz 1 the cross 
section becomes independent of the mediator mass. At 
the same time, for fixed e, the cross section falls as M§,. 
Notice that a similar effect has already been noted in 
dark matter monojet searches [8]. 

For very light mediators, although monojet constraints 
become trivial, other bounds become relevant, for exam- 
ple, rare decays and reactor bounds [ID]. In addition, 
fixed target experiments have been proposed as a probe 
of generic models with light mediators [H]. Finally, star 
cooling bound on NSI may need to be re-examined. 

Our next task is to convert these results into concrete 
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FIG. 6: The acceptance as a function of the Z' mass. Here the 
acceptance is the fraction of events in the initial px > 50 GeV 
parton-level sample that pass the Pythia-level analysis cuts. 
Lighter mediators produce fewer high-pT events, resulting in 
a suppressed acceptance. This is especially evident in the 
veryHighPt case where the choppiness of the curve is a result 
of low statistics. 
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FIG. 7: Our NSI bounds in the Z' model, using CDF and AT- 
LAS monojet data. In the contact limit, the best bounds come 
from the LHCs hardest cuts, while below M z > < 200 GeV 
CDF's softest (GSNP) cut is more constraining. In general, 
the optimal cut is a function of the mediator mass. Recent 
CMS monojet data at 4.7 fb _1 [31] provide a constraint very 
similar to the veryHighPt ATLAS curve shown here. 



bounds on the NSI parameters, as was done for contact 
interactions earlier. This means converting the parton- 
level cross sections into simulated jets and applying the 
experimental cuts. Naively, one might model this step 
by a constant acceptance factor, extracted from the con- 
tact operator analysis. In this way, one would obtain the 
bound given by the red dashed contour in Fig. [5] Yet, 
this would be inaccurate, as we find that the acceptance 
is a strong function of the Z' mass. Passing all of our 
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MadGraph points through Pythia, we find that, depend- 
ing on the LHC cut, the acceptance at lower masses can 
be more than an order of magnitude smaller than in the 
contact limit. 

This is illustrated in Fig. |6j where we plot the accep- 
tance of the three LHC cuts (see Sec. Ill A ) as a function 
of Mz>i fixing e = 0.1. We also show the acceptance 
curves for the two Tevatron cuts. These were obtained 
in a similar way: simulating pp collisions with MadGraph 
on a log-spaced grid of Mz> and e values, passing the 
results to Pythia and finally applying the cuts. 

We present the final bounds in Fig. [7J The results 
are very instructive. While at high masses, Mz' ^ 200 
GeV, the HighPT and veryHighPT cuts from the ATLAS 
analysis give the best bounds, at lower masses these cuts 
become less optimal than the LowPT cut. What is more, 
the best bound in this case comes from the CDF GSNP 
cut, the softer of the two Tevatron cuts. This finding is 
consistent with the decision by the CDF collaboration to 
use the same GSNP cut in measuring the invisible width 
of the Z 

These results demonstrate that for each value of Mz 1 
there is an optimal set of cuts for CDF and ATLAS. 
Then, to achieve maximal sensitivity throughout the en- 
tire Mz> range, both collaborations should vary these 
cuts as a function of the values of Mz> ■ 

Comparing the parton-level results in Fig. [5] with the 
bounds in Fig. [7] we notice that the main effect of includ- 
ing the mass-dependent acceptance is to further weaken 
the sensitivity at low masses. The bound for Mz' ^ 30 
GeV is weaker than in the contact limit. (Low-mass me- 
diators have a more sharply falling px spectrum at the 
LHC and Tevatron compared to intermediate and high 
mass mediators, resulting in lower acceptance.) This pro- 
vides another important reason to go beyond the contact 
interaction limit: neutrino NSI could be mediated by a 
relatively light sector (see also [35])- 

It should be noted that the possibility of such a light 
sector has recently sparked a great deal of excitement, 
in the context of "nonstandard" dark matter physics j£, 
EU 144-49 . We find the possible connections between 
neutrino oscillation and dark matter anomalies intriguing 
and suggest that they be explored further. 



V. DISTINGUISHING NSI FROM DARK 
MATTER 

We have seen that neutrino NSI and dark matter pro- 
duction have similar signatures in monojet events. Es- 
pecially for the flavor off-diagonal NSI, which do not in- 
terfere with the SM, our NSI analysis could be directly 
translated into the language of dark matter. In partic- 
ular, the previous Section shows how the acceptance of 
the ATLAS and CDF experiments to dark matter pro- 
duction events varies with the mediator mass and how 
the experimental cuts could be optimized as a function 
of this mass. 



In order to distinguish neutrino NSI from dark matter, 
or other exotic invisible particles, one should go beyond 
monojets. Specifically, one can search for corresponding 
anomalies in processes involving charged leptons, exploit- 
ing the fact that neutrinos are related to charged leptons 
by SU (2) l gauge invariance. Below we consider a couple 
of examples of such signatures. 



A. Leptoquarks as NSI mediators 

The neutrino-charged lepton connection is, in gen- 
eral, model-dependent. Let us consider, as the first 
example, the specific leptoquark model mentioned ear- 
lier (in Sect. iLVl). The mass splitting of the SU(2) L 



leptoquark doublet could suppress the charged lepton 
rates. It, however, contributes to the oblique T pa- 
rameter [5U] and hence can be probed by precision elec- 
troweak data. Ref. [55] derives the constraint on the ratio 
(Mi/M 2 ) 2 < 5.2 and thus concludes that the neutrino 
NSI in this model can be at most a factor of 5.2 greater 
than the corresponding charged-lepton NSI. The flavor 
changing charged-lepton NSI are, in turn, constrained by 
precision rare decay measurements, such as t~ — > e~7r 
or t~ — > e~rj |14) . The resulting bounds on the charged- 



lepton NSI are found in 



to be e CL < 10" 2 . This 



corresponds to the neutrino NSI bound e eT ^$ 0.05, close 
to the level suggested by solar neutrinos. 

In fact, modern data restrict this model further. The 
constraints on the relevant r decay branching ratios have 
been improved by the BELLE collaboration, by a factor 
of 46 on r~ — > e~ir and a factor of 90 on t~ — > eTr\ [51] . 
Moreover, the mass splitting is well constrained by the 
recent LHC data and updated bounds on the T param- 
eter. The current best fit Higgs mass is now close to 
~ 120 GeV, and LHC data exclude leptoquark masses 
below « 650 GeV [52]. (By comparison, Ref. [26] con- 
siders 95 GeV < M H < 1 TeV and assumes the lightest 
leptoquark mass to be just above Mz/2.) Using a refer- 
ence value of rriH,ref — 120 GeV and that S and U are 
tiny in this model, one finds T < 0.16 at the 95% C.L 
[S3], from which we obtain (Mi/M 2 ) 2 < 1.2. As a result, 
the neutrino NSI parameters in this model are now re- 
stricted to the sub-percent levels, below the sensitivity of 
the solar neutrino data. 



B. Multileptons at the LHC 

As our second example, we consider scenarios in which 
the NSI are described by the dimension-8 operators in Eq. 
(|| up to LHC energies, and estimate their contribution 
to the 3- and 4-lepton final state processes [29"] 



pp^ w + w~e+£p. 



(7) 



Here, the unphysical components of the Higgs in Eq. ^ 
have been "eaten" by the longitudinal components of the 



Multilepton CMS [56] 


3-1 


A-t 


e uP 


uP 


e uP 




e uP 


uP 


e uP 


e uP 


0.19 


0.13 


0.2 


0.5 


0.34 


0.24 


0.36 


0.85 



TABLE II: Bounds on contact NSI couplings using the CMS 
multilepton data, based on C — 4.98 fb [55]. Here i,j = 
e, fj,. All bounds correspond to 95% C.L. This Table assumes 
only one coefficient at a time is turned on. 



W's, as is clearly seen in the unitary gauge. The desired 
signals are obtained when at least one of the W's decays 
leptonically. 

To derive bounds, we use the recent CMS search with 
Js = 7 TeV and 4.98 ftT 1 of data [56] . Our results are 



summarized in Table [TT] and the details of the analysis are 
presented below. Comparing Tables |l| and |Ll| we see that 
the multilepton and monojet searches at present yield 
comparable constraints on the contact operator 

It should be noted that this parity does not hold for 
lower mediator masses. In contrast to monojets, the mul- 
tilepton signal is in fact always weakened when the me- 
diator can be kinematically produced, since there is no 
resonant enhancement. Nevertheless, the multilepton fi- 
nal states represent a very distinctive signature of NSI, 
and should be pursued further. 

As a side note, CMS in fact sees a modest excess 
of multilepton events in certain channels. Intriguingly, 
these events can be accounted for by values of NSI close 
to those suggested by the solar data, e.g., e eT ~ 0.08, as 
described below. 



1. Analysis details 

The CMS multilepton search divides its analysis into 
separate categories depending on the number of hadronic 
ts (N(rh)) identified. Since events with N(rh) > have 
a sizable background, we only consider final states having 
N(rh) = 0. Such final states occur not only for a, (3 = 
e, /i, but also if at least one of the flavor indices is a 
t and any primary rs decay leptonically. Such events 
contribute to the 3£ and 4£ signatures. 

The dominant SM backgrounds for both 3£ and 4£ fi- 
nal states come from Z/7*+jets, it, and vector boson 
production [5H [55] . The primary uncertainty is due to 
the simulation of these backgrounds, and is currently sys- 
tematics dominated. The background from heavy meson 
decays is also relevant to the 3£ search [5 7) . 

Multilepton events satisfying the lepton triggers and 
basic object selection have exactly 3 or 4 leptons, where 
the pt of each lepton satisfies a cut that depends on its 
rank. The highest px lepton has px > 20 GeV, the sec- 
ond has pt > 10 GeV and additional leptons have px > 8 
GeV. In events passing the single-lepton trigger and basic 
object selection the highest px lepton satisfies px > 35 



GeV if a muon, or px > 85 GeV if an electron. Ad- 
ditional leptons in the single-lepton trigger also satisfy 
Px > 8 GeV. Finally, all leptons must be central with 
< 2.1. Leptons are also required to be separated, 
with a separation larger than AR = 0.3 from any nearby 
jet. Jets have px > 40 GeV, are central with |?7| < 2.5, 
and are separated from each other with AR > 0.3. Ad- 
ditionally, to maximize the significance of the signal we 
veto events in which the invariant mass of any opposite- 
sign, same-flavor lepton pair is less than 105 GeV. Similar 
cuts have been imposed by ATLAS and CMS in [531155] , 

We use MadGraph_v5 to simulate the signal, then use 
the built-in pipeline to pass events to Pythia for show- 
ering, ISR, and jet clustering. Since additional leptons 
must have a px larger than 8 GeV, to obtain accurate cov- 
erage of the signal we lowered the generator-level charged 
lepton px cut to 0. For the 3£ analysis, we also lowered 
the generator-level ARjj cut to zero, as discussed below. 

For definitcness, we compute the number of multilep- 
ton events when the monojet bound is saturated. Our 
results are normalized to £ , as given in Eq. This 
accounts for the possibility of multiple NSI couplings si- 
multaneously turned on. 

We first consider the 4£ final state, where both W's 
in Q decay leptonically. We find that after the event 
selection described above, A£ signal events are dominantly 
(60%) in the high MET > 50 GeV, low H T < 200 GeV 
region "(high, low)". After all object selection the signal 
cross section in this (high, low) region is found to be 
approximately an = 0.3£ fb, corresponding to 



N 4 



1.5 £ , i, j — e, [i. 



(8) 



We now turn to leptonic rs. Simulating r decay is 
beyond the scope of the present work. Instead we use 
the same efficiency of « 60% (not including the leptonic 
BR) as we found for the ee/e/Lt//z/z-type NSI. This esti- 
mate of the efficiency is reasonable as the primary rs are 
centrally produced and their px spectrum has a median 
at ~ 400 GeV, with over 90% of the rs lying above 100 
GeV. Therefore, a lepton produced from the decay of a 
primary r should have a px value large enough on aver- 
age to pass the multilepton triggers. Moreover, here the 
MET distribution is expected to shift to higher values 
compared to ee/e/x//z/z-type NSI, for which ~ 80% of the 
signal lies above MET = 50 GeV. We find 



( , J 0.5 £, i = e or fj,, 
1 0.2£, i = r. 



(9) 



In the "no Z, (high, low)" region of the U sample, CMS 
observes a single event, with an expectation of 0.2 ±0.07. 
Using Poisson statistics, 5 signal events are allowed at 
95% CL, resulting in the numbers in Table [n] Compared 
to the monojet bounds in Table [TJ the constraints from 
the A£ analysis are weaker by a factor of ~ 2 on the e ee , 
£ efl and e MM couplings, by a factor of ~ 3 on e eT or e 
and by a factor of ~ 5 on s TT . 
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We now turn to the 3-lepton final state, which occurs 
if either one of the Ws in Q decays leptonically. We 
neglect the contribution of primary 4-lepton events, in 
which one of the leptons does not pass the event selec- 
tion. One characteristic of the process ^ is that the Ws 
are boosted, with a median pt — 300 GeV and ~ 80% 
having px > 150 GeV. Due to the boost, the two partons 
produced in the hadronic decay are collimated, with a 
median Ai? ~ 0.6. For optimal coverage of the signal in 
our MadGraph event generation, we therefore lowered the 
generator- level Ai? separation between partonic jets to 0. 
(This is possible because the partons from the decay of a 
boosted W do not suffer from a collinear singularity.) At 
the analysis level we then find the ARjj > 0.3 separation 
cut retains over 99% of the signal. Imposing the selec- 
tion criteria as described above, we find that ~ 60% of 
the signal is in the (high, high) bin. In this bin, we find 
after all object selection a^i = l-8£ fb, corresponding to 

!9£, i,j = e,fx 
3£, i = e,mj = T (10) 
1£, i = j = r. 

CMS observes 8 events in this bin, with an expectation 
of 5±1.3 events. Using Poisson statistics, 11 signal events 
are allowed at 95% CL. Proceeding as before, we find 
that Eij < 0.19 for diagonal NSI ee or fifi, and e eM < 
0.13. Here the bounds from the 3£ final state are only 
slightly weaker than those from the mono jet analysis. 
For processes producing a primary r we make the same 
approximations as in our At analysis, and only consider 
the contributions from leptonic rs. Then e e /^ T < 0.2 and 
e TT < 0.5 at 95% CL. 

One can also view the multilepton data in a different 
light and ask to what extent they prefer nonzero NSI 
couplings. Intriguingly, the NSI couplings favored by so- 
lar data may provide a better fit to the multilepton data 
in 56 . For example, if e er = 0.08, NSI will contribute 
~ 1 event in the high MET, high bin for the 3 lepton 
case, but only ~ 0.2 events in the high MET, low Ht 
bin for 4 leptons. CMS sees excesses of 3 and 1 events 
in these two categories. More statistics will be required 
to determine whether these excesses are truly due to new 
physics or simply upward fluctuations. 

Multilepton events at the 8 TeV LHC run will probe 
the contact operator in Eq. ^ even further, since the in- 
crease of energy from 7 TeV will double the signal cross 
sections. With the same object selection criteria, we find 
a u = 0.72£ fb and a u = 4.5£ fb at 8 TeV. Normaliz- 
ing the NSI couplings e,-j (for i,j = e, fx) to the mono jet 
bounds and using the same selection criteria described 
above, then with an integrated luminosity of 20 fb -1 
there will be ~ 16 4-lepton and ~ 90 3-lepton signal 
events in the (high, low) and (high, high) regions respec- 
tively. Alternatively, assuming e eT — 0.08, these regions 
will have ~ 2 4-lepton signal events and ~ 14 3-lepton 
signal events. 



VI. CONCLUSIONS 

In this Letter we proposed using the monojet plus miss- 
ing energy datasets at the Tevatron and the LHC as a 
novel probe of nonstandard neutrino interactions. As- 
suming first that the NSI remain contact at the LHC 
energies, we derived stringent bounds on the parameters 
egf, e?£, and e%l with q = u, d and P = L,R. These 
bounds come from ATLAS's 1 fb -1 dataset [31], which 
has already overtaken the Tevatron's CDF experiment in 
sensitivity in this regime. The bounds are summarized in 
Table [I] and approach (within a factor of 2-4) the levels 
motivated by the solar neutrino data. 

Given this state of affairs, further progress is highly de- 
sirable. We note, however, that the present bounds, while 
based only on only on 1 fb -1 of data, are already system- 
atics dominated. Further improvement in NSI monojet 
bounds is therefore largely predicated on improving our 
understanding of the systematics at the LHC. 

Our monojet bounds apply to neutrino-quark interac- 
tions in a flavor-independent way since the processes in 
Fig. [2] are neutrino flavor-blind. Importantly, they also 
apply equally well to sterile neutrinos with couplings to 
SM quarks |46j . or to light dark matter models as dis- 
cussed in Sect. |V] 

We also considered the effect of relaxing the contact op- 
erator assumption, thereby allowing the mediator of new 
physics to be directly accessible at current LHC energies. 
In this case, the analysis inevitably becomes model de- 
pendent. We showed that with an s-channel mediator 
the bounds are particularly stringent if the scale of new 
physics is in the range of ~ 10 2 GeV. At the same time, 
new physics below < 30 GeV could escape the mono- 
jet bounds and appear first in neutrino oscillation ex- 
periments. Thus, NSI with observable oscillation effects 
could originate either at high scales, > 2 TeV, or in the 
low mass window, < 30 GeV. 

It is noteworthy that light mediators have recently at- 
tracted considerable attention in connection with models 
of dark matter gT] I41H55] . We find interesting the 
possibility that the solar neutrino data may also favor 
new physics at the same scale. 

We have seen that in the regime < 200 GeV CDF with 
its soft GSNP cuts actually bests ATLAS in its NSI sen- 
sitivity. We encourage the Tevatron and LHC collabora- 
tions to publish their monojet results with an extra low 
cut, or perform the analysis of neutrino NSI themselves 
tuning the cuts as a function of the mediator mass. 

The finite-mass mediator scenario also allows us to ad- 
dress a posteriori the range of applicability of the earlier 
contact interaction analysis. As is evident from Fig. [5] 
the contact limit sets in only for masses above several 
TeV. Physically, this means mediators lighter than that 
may be produced directly in high-energy collisions. We 
stress that the analysis in this case is by necessity model- 
dependent. This scale will be pushed up even higher as 
the energy of the LHC beams is increased. 

The present monojet dataset provides bounds si- 
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multaneously on neutrinos, dark matter, and extra- 
dimensional models. While many of the analysis steps 
are similar, there are several important distinctions of 
neutrino NSI compared to the other two types of new 
physics. First, flavor-diagonal NSI interfere with the SM 
processes. Therefore, further experimental improvements 
can lead to much more stringent limits (linear in e's). 

The second important difference is that neutrinos are 
part of an SU(2)l doublet and hence can contribute to 
processes involving charged leptons. Signatures in the 
monojet and multilepton search channels are thus cor- 
related. We have considered an example of this in 



Sec. VB where we find the bounds from 3-£ multilep- 
ton events on NSI couplings of the first and second gen- 
erations are practically identical to those from mono- 
jets. Using the monojet bounds derived in the earlier 
sections, and the values of the NSI parameters hinted 
at by the present-day solar neutrino data, we found a 
predicted multilepton rate that is curiously close to the 
just-reported hints of excess [55 . Multilepton searches 
at the 8 TeV run of the LHC will probe NSI even further. 

Lastly, for the finite mass scenarios, the best way 
to search for the physics behind neutrino NSI becomes 
model-dependent. The monojet analyses should then be 
viewed as part of the NSI search portfolio, providing di- 
rect though not necessarily strongest bounds. Even from 
the limited discussion here it is clear that such seem- 



ingly disparate searches for leptoquarks, Z''s, multilep- 
tons, and monojets could have a connection to each other 
and to the data in neutrino oscillation experiments. We 
urge the LHC collaborations to seriously consider a co- 
herent program targeting neutrino NSI physics with mul- 
tiple search modes. We ourselves plan to return to this 
problem in a future work. 
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Note added: The present updated version of this Letter 
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neutrino survival probability in Fig. [T] with the measured 
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the Higgs searches in our discussion of the leptoquark 
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